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Abstract. Transient spectral features have been discovered in the X-ray spectra of Active Galactic Nuclei, mostly in the 5-7 keV 
energy range. Several interpretations were proposed for the origin of these features. We examined a model of Doppler boosted 
blue horns of the iron line originating from a spot in a black hole accretion disc, taking into account different approximations 
of general relativistic light rays and the resulting shift of energy of photons. We provide a practical formula for the blue horn 
energy of an intrinsically narrow line and assess its accuracy by comparing the approximation against an exact value, predicted 
under the assumption of a planar accretion disc. The most accurate approximation provides excellent agreement with the spot 
orbital radius down to the marginally stable orbit of a non-rotating black hole. 
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1. Introduction 

X-rays are the main source of information about the inner- 
most parts of accretion discs in active galactic nuclei (AGN) 
and Galactic black hole systems. The iron Ka line is the most 
prominent spectral feature, occuii'ing at a rest energy Eq rang- 
ing from ~ 6.4 keV (neutral iron) to ~ 6.96 keV (H-like iron). 
The natural width of these lines is only a few eV (Thompson et 
al.lMm. They are thought to originate as a result of illumina- 
tion of an accretion disc, which is formed by relatively cold ma- 
terial, by high-energy continuum photons from a corona or a jet 
above the disc (Fabian et al. I2000I Reynolds & Nowak I2()()3> . 
The Ka lines are multiplets. The neutral iron line is a doublet 
with energies Fkq-i = 6.404 keV, Fkijj - 6.391 keV. The small 
energy difference makes the doublet unresolved with present 
spectrometers, though it may be resolvable by ASTRO-E2 in 
the near future. Another prominent spectral feature, namely 
the K/3 line, arises at an energy of ~ 7.07 keV and the iron 
edge occurs at 7.1 keV. Special and general relativity effects 
strongly modify the appearance of all these features. For in- 
stance, the continuum becomes broader with increasing disc 
inclination due to enhanced Doppler shifts. The photoelectric 
edge is blurred into broad troughs. X-ray spectroscopy can then 
be employed to reveal the effects of general relativity influenc- 
ing the observed energy and thereby shaping the spectral lines 
and the continuum slope. 

Fast circulation of gas is the dominant source of line broad- 
ening and forms, together with substantial gravitational red- 
shift, the observed profile. A double-horn shape is expected 
in the simplest case of a geometrically thin, Keplerian disc 



(Fabian et al. I1989t . and the discovery of this kind of profile 
has been indeed reported in several objects, the best-known 
examples being the Seyfert galaxy MCG-6-30-15 (Tanaka et 
al. 1 1995b and the microquasar GRS 1915-1-105 (Martocchia 
et al. lTOOl . where the line shows equivalent widths reaching 
~ 200 eV. It has been argued that fine substructures of this 
broad feature could be used to constrain the inclination angle 
of the source, radial emissivity in the disc plane, and even the 
angular momentum of the central black hole (see Martocchia et 
al. 120001 Beckwith & Done l2004l and references cited therein). 
However, very broad “classical” double-horn relativistic lines 
are often absent in AGN (e.g. Bianchi et al. 12004b . and it is 
important to seek additional evidence. 

In recent years, narrow emission features have been re- 
ported in the X-ray spectra of several AGNs: see Turner et al. 
( l200l l2004l l2005t : Guainazzi (|2003; Yaqoob et al. (|2003l; 
Bianchi et al. J2004b : Dovciak et al. J2004ab : Gallo et al. 
J2004b : McKernan & Yaqoob J2004b : Porquet et al. J2004b : 
Della Ceca et al. J2005b . This topical discovery renewed the 
interest in accretion disc hot spots. The narrow lines we have 
on mind are typically found in the 5-7 keV energy range, i.e. 
mostly on the redshifted side of the iron line rest energy. A ten- 
tative explanation for these puzzling features, which we adopt 
in this paper, is in terms of iron emission, which originates in 
a localized area of the disc (6r < r^). According to this scheme 
the spectral features are produced by transient magnetic flares 
briefly illuminating the disc and producing the iron line by flu- 
orescence in the surface layer (Haardt et al. 1 19941 Poutanen & 
Fabian[l992 Czerny et al. 12004b . 
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The idea of revealing the signatures of black holes via spec- 
troscopy of orbiting spots has great potential and a long history 
(Cunningham & Bardeen [l973> . Theoretical line profiles have 
been studied including specific effects of general relativity: 
the frame-dragging (Koiima [l991l Laor ll99Tl Bromley et al. 
nw7t . extreme light-bending and multiple images (Matt et al. 
fT9^ Viergutz[T993j Bao et al.[T994l Zakharov & Renin lMBl 
Beckwith & Done 1200^ . non-axisymmetric (spiral) waves in 
accretion discs (Karas et al. l2001l Hartnoll & Blackman l2002t . 
and self-gravity (Karas et al. Usui et al.[T99l. Various 
applications have been proposed for the general-relativity en- 
ergy shift of spectral features from co-rotating spots; in particu- 
lar, it was suggested that the range of energy fluctuations could 
constrain the source geometry and discrimate between planar 
and thick discs (Bao & Stuchlik [l993> . 



In this paper we examine extremal values (maximum and 
minimum) of observed energy of a line from an orbiting spot. 
The main purpose of this study is to provide a simple, practical 
and accurate formula that will be useful in analyzing the red- 
shifted narrow lines. We neglect light bending, because its in- 
fluence is important for only a relatively minor fraction of pho- 
tons on special rays, i.e. those passing near caustics. This ap- 
proximation was assessed by Gerbal & Pelat J199H . but these 
authors did not aim to provide an explicit formula that would 
be convenient for practical purposes. Our approach is similar 
to Zhang & Bao J199U , who also ignored light bending. We 
improve their treatment of light aberration and this allows us 
to express the redshift function (defined as a ratio of observed 
to emitted energy; see eq. i} below) in a way that is compara- 
bly simple to but more accurate than their eq. (2.9). As a fur- 
ther step, we include light bending in terms of Beloborodov’s 
<2002t approximation and reach fractional accuracies of a few 
percent or better for the redshift function; the extremal values 
are given with relative error better than 0.1%. Thus the result- 
ing shift of energy is estimated at least an order of magnitude 
more precisely than before, and the error of its extremal values 
is negligible in most practical cases. 



In the next section we summarize the main results of this 
paper. In particular, in subsec tion l2.1l we specify our assump- 
tions and give a prescription for calculating the observed en- 
ergy range spanned by the spectral line (i.e. extremal values of 
g-factor). It will be of interest to a reader seeking a straight- 
forward approach to practical problems of data fitting. Then, in 
subsec. o we describe the individual steps of the derivation, 
including a comparison with previous results of other authors. 
In Sec. 13 in a brief re-discussion of several AGNs we illus- 
trate the interpretation of narrow spectral features in terms of 
the blue horn produced by an orbiting spot. The results pre- 
sented here are limited to static black holes and geometrically 
thin, Keplerian discs. Even with these limitations, they may be 
useful in interpreting the present data and deriving disc param- 
eters. As it will be shown in section|3 the narrow features de- 
tected so far do not require spinning black holes. 



2. Calculation of the redshift 

2.1. Summary 

For the gravitational field of the system we assume a non- 
rotating black hole and so adopt the Schwarzschild spacetime 
metric (Misner et al. ll972> . 

di^ = df^ H dr^ + (d6^ + sin^fl d(^^) , (1) 

r r -2 '' ' 

where the radius r is evaluated in units of gravitational radius 
Tg s GMc~^^\.5 X 10^(M/Mq) cm, with M being the mass 
of the black hole. In eq. o as well as everywhere in this text, 
geometrized units with c = G = 1 are used; distances are thus 
made dimensionless by scaling them with M. Notice that the 
accretion disc own mass can be neglected, i.e. M^<s:M, because 
we constrain ourselves to the innermost regions of the flow 
where the disc is non-self-gravitating (in terms of Toomre’s cri- 
terion). 

The disc is assumed to be Keplerian and extending down 
to the marginally stable orbit, r^s = 6, below which it ceases 
to exist. An observer is located at infinity on the x axis (i.e. 
(j) - 0°). By approximating trajectories of light according to 
Beloborodov J2002t and evaluating the effect of aberration as 
detailed in subsec. 12.21 we obtain the extremal values of the 
energy shift in the form 

r -ni/2 ■ ^ 

r + r-2-i-4(l + cos ([>± sin 0 q) sin ([>± sin 

These extremal values occur for the azimuthal polar coordinate 



cp = (p±, defined by 

/ 7 \l/2 

sin0± = 4- ( 1 - cos , (3) 

COS,^* = /r‘ +(r- - 2(r - 1)] (4) 

(<p = 90° corresponds to the maximally receding part of the 
disc). The functions fi(r, 6o) and / 2 (r, 6o) are 

fi(r,9o) = 3(r-2)sin6>o, (5) 

fiir, 0o) = (r - Af + / 3 ^ + [2 (r - 4)^ + /a, (6) 

with /a(r, 6»o) being 

Mr, 6o) = 3^ (r- 2) cos 0o- (7) 



Equation d3 can be used to determine the range of energy 
spanned by a line without the need to evaluate special functions 
(elliptic integrals and their inverse in case of Schwarzschild 
black hole; e.g. Chandrasekhar IT^^ Cadez et al. l2003> or to 
perform numerical ray-tracing in the curved spacetime. Instead, 
purely algebraic formulae are involved, allowing us to give ex- 
tremal values of redshift as functions of the spot radius r and 
observer’s inclination angle Oo- Below we estimate the accu- 
racy of our approximation and And it very precise even at large 
values of inclination (almost edge-on, i.e. 0^ ~ 89°) and small 
radii (down to the marginally stable orbit, r ~ r„is). We And 
much better agreement between approximate and exact values 
than what can be achieved with the Zhang & Bao (I1991t ap- 
proach. 
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Fig. 1. Accuracy assessment for different approximation formulae, (i) Top row - left panel; the dependence of the maximum and 
minimum relative error 6±g{r, Og) of the .g-factor on observer’s inclination angle Og for a narrow ring, r = r„is- Solid lines refer 
to eq. (O (the upper one is for maximum, the lower one is for minimum); the dotted line refers to eq. d; the dashed line is 
for eq. (O- The approximation formula J1 2\ is the worst of the three cases in the sense that it has a signihcant non-zero value 
of the minimum error. This is due to neglecting aberration. The relative error 6g - Igappiox - g\!g is determined with respect to 
exact (numerically computed) values, (ii) Top row - right panel: the relative error 6g±{r, 9g) of the extremal values of g-factor for 
a ring r - rms as a function of observer’s inclination angle. The relative error dg± = (^approx “ g±)lg± of the extremal values is 
again determined with respect to exact values from equation of geodesic. The solid/dashed lines are for maximum/minimum of 
g-factor; both lines are for the approximation O- (iii) Bottom row - left panel; the same as top right but with the approximation 
O- (iv) Bottom row - right panel; the same as top right but for the approximation d3. 



2.2. Details of the calculation 



In this section we derive the results and provide comparison 
with previous works. The relativistic shift g of the energy of 
photons, i.e. the ratio of the frequency observed by a distant 
observer to the frequency emitted in a local frame co-rotating 
with the disc material, can be expressed, within the approxima- 
tion of geometrical optics, in terms of the four-momentum of 
photons p and the four-velocity U of the emitting material and 
of the distant observer. The shift g(r, (p\ 9g) is a function of polar 
coordinates in the disc plane and of observer’s inclination. 



gfr, 4>; 9g) 



Eg 



PofjUo 

PgaU^ 



Pot 

PgaU^ ’ 



( 8 ) 



where Eg denotes the photon energy measured by a distant 
static observer (r ^ oo). Eg is the locally emitted energy from 
a source co-rotating with the disc, Ug is the four- velocity of the 
disc medium, and Ug is the observer’s four-velocity. 

Equation (|H|l can be written in the following way; 

g(r, 0; 0o) = f T , (9) 

U‘g(r)[l-D(r)l(r,tf,-9g)] 

where t/'(r) = (1 - 3/r)“'^^ is a component of the disc 
Keplerian four-velocity, Q(r) = is the angular veloc- 

ity and l(r, cj)', 9g) is the azimuthal component of the photon’s 
angular momentum (constant along null geodesics). The de- 
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Fig. 2. Left; the dependence of the maximum and minimum relative error S±g(r, ffo) of ^-factor on the radius for inclination angle 
ffo = 60°, in logarithmic scale. Solid lines are for the approximation O (upper for the maximum, lower for the minimum), 
dotted line is for the approximation J14> and dashed line is for the approximation O- Right; the relative error Sg±(r, Oq) of the 
extremal values of g-factor. The solid lines are for eq. (inii and the dashed lines are for eq. lll5> . Upper lines show the relative 
error of g+, lower lines are for the relative error of g-. The relative error for the approximation 0 is better than 0.1% and thus it 
is not plotted here. 





Fig. 3. Left; g-factor as a function of the azimuthal coordinate (p for r - rms, 9o - 60°. The solid line is for eq. C3, the dotted 
line is for eq. (O and the dashed line is for eq. O- The exact numerical solution is also plotted by the dashed line and cannot 
be distinguished from the approximation dn} in this figure. Right; as on the left but for the relative error of 6g for the three cases 
of approximation formulae. Notice that the error of the best one, eq. o. is virtually zero. 



pendence l(r, cf>', 6o) on the polar coordinates in the disc can 
be calculated analytically with the help of elliptic functions 
( Chandrasekhar fT992> . Let us adopt the approximation in terms 
of the 0-component of the photon momentum in a local frame. 



l(r, 0; 6»o) 



0; 0 ^) -I- 

nP‘f'\r, 0; 6»o) 4- 






where is a unit three-vector in the direction of 

photon emission. We choose the local frame attached to an ob- 
server with the tetrad e" ^ (here, we assumed that p, - and 
p<p = 0- 

The simplest assumption one can impose is that photons are 
emitted “straight” toward the observer, 

n^^\r, 0; do) - - sin 0 sin Oo- 



( 10 ) 



(11) 
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Fig. 4. Left: the maximum redshift factor g+(r) as a function of the radius. Different curves are parametrized by inclination angle 
Oo- Right: the radius at which the function g+(r) is maximum as a function of Oo- For 60° the maximum occurs at marginally 
stable orbit where the inner edge of the disc is assumed. 



This was adopted by Zhang & Bao IMB and in many subse- 
quent papers, neglecting the effect of light bending. This turns 
out to provide a sufficiently accurate description in many astro- 
physically relatistic situations; however, Zhang & Bao further 
assume that the photons move straight in the direction to an ob- 
server in the local frame co-moving with the disc (i.e. the tetrad 
in eq. (iin) corresponds to a Keplerian observer and vector 
is expressed in this tetrad). In this way they also neglect 
the aberration due to rotation of the disc, which enlarges the er- 
ror especially at small inclination (pole-on view). The resulting 
formula for the redshift is 



gzBir, (jy, Oo) 



r-2 

r(r-3) 



1/2 

|^(r - 2)*^^ - sin0 sin0o] ■ 



( 12 ) 



Extremal values occur on a ring r - const (p = - -i-90°. 



^zb(t Oo) 



r-2 

r{r-3) 



1/2 

[(r-2)*^^ + sin6»o] . 



(13) 



neglecting the light bending. In this approximation the tetrad of 
a static observer is used in eq. (0 and the redshift function is 



gsi(L 0; Oo) = 



(r-2)(r-3) 



1/2 



(r - 2)F2 + sin (^sin0o ’ 



(14) 



and its extremal values on a ring, which again occur at (p± - 
+90°, are 



^fi(T ^*0) = 



f(r-2)(r-3) 



1/2 



1 

(r - 2)F2 + sin0o 



(15) 



This approach is very successful at estimating the extremal val- 
ues (see the left panel in bottom row of Figure^ but the relative 
errors of the redshift factor itself remain non-negligible. 

Finally, the approach can be refined by assuming that pho- 
ton trajectories are curved by gravity and adopting a suitable 
approximation to their shapes (Beloborodov ’s|2002|l. Then the 
azimuthal component of the directional three-vector is (see 
AnnendixiAl 



These approximate formulae suffer from excessive inaccuracy 
in some regions of the parameter space. We computed these 
fractional errors and illustrate them in three plots. In Figure ^ 
we plot the maximum and minimum relative errors of the red- 
shift factor 6±g(r, Oo) and the relative error of the extremal val- 
ues of g-factor 6g±{r, Oo) on a ring r = rms (marginally stable 
orbit). Figure 12 shows the dependence on radius for Oo = 60°. 
Figure|2shows the azimuthal dependence of g-factor itself and 
its relative error. 

The main source of error is not the neglect of light bending 
and, indeed, considerably better results are achieved if we per- 
form this calculation in the local frame of a static observer, still 



n^‘^\r, cf>-, Oo) - — 



1 - cos^ a 

1 - cos2 l/f 



1/2 

sin0sin0o, 



(16) 



with cos q; = l-(l-cosi/f)(l-2/r)andcosi/r = cos0sin0o- The 
angle a is the angle that local static observer measures between 
the radial direction and the direction of emission. The redshift 
factor is 



gB ( r ,( p ; 0o ) 



[r(r-3)] 



1/2 



r _i 1 1/2 

r + r - 2 + 4 (1 + cos (^sin0o) sin(^sin0o 



( 17 ) 
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Source 


Elme [keV] 


r [GM/c^] 


[deg] 


^0o=const [C?A//c“] 


Reference 


NGC 3516 


5.57 


6-12 


0-23 


6.6 (Bo = 20°) 


Turner et al. i2002t 




6.22 


6-50 


0-35 


12.4 (00 = 20°) 






6.53 


>6 


20-40 


56/93 (00 = 20°) 




ESO 198-G024 


5.70 


6-14 


0-26 


7.2 (0o = 20°) 


Guainazzi <2003 > 




5.96 


6-21 


0-30 


8.9 (0o = 20°) 


Dovciak et al. <2004al 


NGC 7314 


5.84 


6-17 


0-28 


6.3 (00 = 27°) 


Yaqoob et al. ^20031 




6.61 


>6 


27^1 


20/93 (00 = 27°) 




Mrk 766 


5.60 


6-12.5 


0-24 


6.7 (0o = 20°) 


Turner et al. j2004J 




5.75 


6-15 


0-27 


7.4 (0o = 20°) 




ESO 113-GOlO 


5.40 


6-10 


0-20 


6 (00 = 20°) 


Porquet et al. <2004> 



Table 1. The narrow features detected so far in AGN which can be interpreted as the blue horns of a 6.4 keV iron line arising 
from an orbiting spot. 



Expressions for the extremal values in this case are summarized 
in the previous section, see eqs. (EJi-Q. 

3. Comparison with observations 

Chandra and 'KMM.-Newton have discovered narrow and 
sometimes transient features in the 5 to 7 keV energy range. 
Dovciak et al. (I2004at argued that these lines may be the blue 
horn of an iron line profile from an emitting annulus in the ac- 
cretion disc, a situation which could arise if an X-ray flare just 
above the disc survives a non-negligible fraction of the orbital 
period or more. 

These features offer the possibility to measure the black 
hole mass (see discussion in Dovciak et al. I2004at . and 
indeed Iwasawa et al. applied this method to the 

Seyfert 1 galaxy, NGC 3516. The quality of data, however, 
is in most cases insufficient for tfiis purpose. For the time be- 
ing, and waiting for the next generation of X-ray satellites like 
Constellation-X and XEUS, we estimate the orbiting spot pa- 
rameters. 

In Figure0](left panel) we show the maximum energy shift 
g+{r), which cotTesponds to the blue edge of the line and the 
peak flux when lensing is neglected or unimportant. Only di- 
rect image photons were considered (higher-order images can 
be important only for equatorial inclinations of 6o ~ 90°; see 
below). In this case the energy shift is given by eq. 0. We 
checked that the exact value (based on the full geodesic equa- 
tion and the ray tracing method) for the energy shift is indistin- 
guishable from the approximation adopted in this figure. This 
plot permits determination of ranges of allowed (r, flo) pairs 
from the measured energy of the narrow feature. A comple- 
mentary plot is given in the right panel where we show radius 
at which the maximum energy shift occurs for a given value of 
inclination angle 9 - 6q. 

One can ask whether the lensing effect, which we have ne- 
glected in our present discussion, could change the conclusion 
about the observed energy of the Doppler boosted horn of a 
relativistic fine. The reason is that for higher inclination an- 
gles (6o ^ 60°) the lensing effect and bending of light rays 
(that changes the projected area of the spot) can play an im- 



portant role: thus the dominating feature need not be a blue 
horn associated with the maximum shift of energy but, instead, 
a new peak of the maximum observed flux that should arise 
redwards of g+(r) (see e.g. Matt et al. ll993> . We have there- 
fore employed a ray-tracing code (Dovciak et al. l2004bl taking 
all these general-relativity effects into account, and examined 
time-dependent model spectra. Indeed a new peak is formed by 
photons arriving from another azimuth (different from -90°) of 
the spot orbit. The secondary maximum is in the point where 
the effect of energy shift, lensing, projection area and time de- 
lay combine to produce the maximum amplification. However, 
the time delay, and especially the function g{r, (p), vary rapidly 
across the region of the disc most affected by lensing. As a 
result, the new peak becomes spread over a broad range of en- 
ergy and does not form a narrow line. In other words, even if 
the lensed feature may have a large equivalent width and can 
contribute significantly to the total photometric flux, this effect 
is not relevant for our present discussion of the narrow line ori- 
gin, as the resulting feature is not natTow at all. 

In Table^we summarize the values of r and 9o for the most 
convincing observational cases known. For all features a solu- 
tion in which the orbiting spot is in the innermost stable orbit 
for a Schwarzschild black hole, i.e r = 6, is acceptable. The 
inclination angle is generally expected to be low or moderate 
(the inner accretion disc is presumably aligned with the outer 
torus). If the features we discuss originate at or near the in- 
nermost stable orbit, the corresponding inclination angle is the 
upper value in the table. Naturally, as this angle must be the 
same for different features of the same source, the actual al- 
lowed range is the intersection of the various ranges: 20-23° 
(NGC 3516), 0-26° (ESO 198-G024), 27-28° (NGC 7314), 0- 
24° (Mrk 766). Assuming, therefore, the same inclination an- 
gle for a given source, and assuming this angle to be 20° for all 
sources but one, the cotTesponding radii are fisted in the fourth 
column of the table. For NGC 73 14, if the two features are both 
the blue horn of an annular emission, then the inclination angle 
is constrained to be between 27° and 28°; for this source an an- 
gle of 27° has been assumed. For the only two features that are 
blueward of the line rest energy, the resulting radii could not be 
determined uniquely. This is explained by the fact that the fine 
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Fig. 5. A schematic representation of the plane of photon motion (left), and the definition of angles and directions in 
Schwarzschild coordinates (r, 9, (p) (right). The direction vector n is tangent to the photon ray projected into three-space. 



of g = const can intersect the curve g+{r) at two different radii 
(see Fig.0}. 

The values of the parameters in the table have been derived 
assuming a rest-frame energy of the line Eq ~ 6.4 keV, as ap- 
propriate for neutral matter. However, the iron may be ionized 
and the corresponding Eq higher. The maximum possible en- 
ergy is for H-like iron, and is Eq ~ 6.96 keV (the line is a dou- 
blet again). Because in this case the redshift has to be larger, the 
upper limits to the radius and the inclination angle are reduced. 
For instance, for the reddest feature (that of ESO 113-GOlO) 
the allowed ranges are now 6 < r < 7.6 for the dimensionless 
radius and 9o < 11° for the inclination. 

As discussed above, for high inclination angles the max- 
imum flux could occur not in correspondence with the blue 
peak, as we assumed in this paper, but instead with the mo- 
ment of the spot crossing a caustic at an intermediate phase. 
In this case the resulting inclination comes out larger. We have 
explored and rejected this lensing effect as the origin of narrow 
lines not only due to the large inclination, but also because the 
width of the resulting spectral feature exceeds substantially the 
width of the blue peak. In no case do the observed features re- 
quire a spinning black hole, as a solution with r > 6 is always 
found. 

4. Discussion 

Variability properties contain a wealth of information on ac- 
creting black holes, but much of it is lost in the process of time 
integration in the detector. Several ways have been proposed 
to measure the parameters of accreting black holes, which rely 
on the source variability: for example, limits on the black-hole 
mass and its rate of rotation can be derived from time delays 
between variations of the emission-line strength and of the con- 
tinuum, and from temporal changes of the observed emission 
lines. Time-dependent spectra have not yet provided the an- 
ticipated breakthrough in our understanding; they still remain 
a challenge for theorists and observers alike. This is caused 
by difficulties and non-uniquenesses in the interpretation of the 
available dynamical spectra. On the side of observation a trade- 



off must be accepted between energy and temporal resolution 
which also compromises the usefulness of timing spectroscopy. 
That is why we search for sharp features, and the narrow red- 
shifted lines are promising. 

Here we showed that a simple and accurate formula can 
be given for the estimate of the energy shift. The analytical 
formula circumvents the necessity of numerical integration of 
photon paths and can be straightforwardly used to determine 
the orbiting spot parameters from the measured energy of the 
narrow feature, with an accuracy that is more than sufficient. 
The use of this formula implies the assumptions that the feature 
corresponds to the blue peak in profile, that the effect of black- 
hole rotation is negligible and that the disc is essentially planar. 
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Appendix A: Direction of photon emission 

In the Schwarzschild metric a ray can be defined by three 
points: the origin of the coordinate system (which coincides 
with the center of the black hole), location of the spot (in the 
disc plane), and the observer (at radial infinity). See Figure 0 
for the geometrical setup. Then the directional vector of a 
photon can be expressed with respect to the local static frame. 
We have employed Beloborodov’s (I2002t approximation of 
light rays. One can write 

- cosae^^j - (1 - cos^o;)'^^e^^j, (A.l) 

where the unit vector is perpendicular to the tetrad vector 
and lies in the plane of the photon’s motion. Thus (see the 
right panel in Figure]^ 

= r"'(l -cos^i/')“'^^(0,0, cos0o,sin0sin0o) (A. 2) 
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with cos Ip - cos^sinflo- The components of the directional 
vector = n^e[a)ii in a local static frame are 



= cos a, 



n 



(9) _ 



1 



• cosher 



n 



( 0 ) _ 



1 - COS^ljj 

1 - cosher 
1 - cos^l/r 



1/2 



cos 6 * 0 , 



1/2 



sin sin 00 V 



(A.3) 

(A.4) 

(A.5) 



Here, we have used the tetrad of a static observer located in the 



equatorial plane: 

= (l-2/r)-i/2(i, 0,0,0), (A. 6 ) 

= (l-2/r)'/2(Q^ 1,0,0), (A.7) 

= r-i(0,0,l,0), (A. 8 ) 

= r-'(0,0,0,l). (A.9) 
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